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SUMMARY 

A method is presented that enables convenient determination of 
the pressure drop sustained by air flowing at high subsonic speeds in 
a constant-area passage under the simultaneous influence of friction 
and heat addition. The method is applicable for arbitrary heat-flux 
distribution along the flow-passage length. Air-pressure-drop working 
charts based on the method developed are presented. 


INTRODUCTION 

Many problems arise in aircraft heat- exchanger practice that 
involve the transfer of heat at high flux rates to a compressible 
fluid flowing at high speeds through constant-area flow passages. For 
example, in some recent heat- exchanger designs the temperature differ- 
ential between heat- exchanger wall and fluid approached values of the 
order of 1200° R and the fluid velocities approached the choke condi- 
tion (Mach number of 1). 

Methods have been developed for determining the pressure drop 
experienced in the passage of air through the radiator tubes of recip- 
rocating engines (references 1 and 2). These methods contain simpli- 
fying assumptions that, although adequate over the range of conditions 
encountered in aircraft radiators, cause appreciable error at the high 
airspeeds and heat- input rates considered herein. 

In reference 3, the basic differential flow equation describing 
the pressure variations of a compressible fluid under the simultaneous 
action of friction and heat addition is numerically integrated for the 
specific case of air heated at constant passage-wall-temperature con- 
ditions. In reference 4, a closed-form solution of the basic differ- 
ential flow equation is obtained for the special case of an exponential 
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variation of fluid temperature with distance along the flow passage. 
Many cases arise, however, where the heat input to the fluid is a 
complicated function of distance along the flow passage. 

In an investigation made at the NACA Lewis laboratory and 
reported herein, the basic differential flow equation is formally 
integrated by a method based on a simplification that is shown to 
introduce little error in the calculation of pressure drop over the 
Mach number and heat-input range encountered in high-flux -rate heat 
exchangers. The integration results are presented for air (constant 
ratio of the specific heats equal to 1.40) in convenient chart form. 
The method is applicable for arbitrary heat-flux distribution along 
the flow passage and for Mach numbers up to approximately 0.90. 

Detailed examples are given that illustrate the methods presented 
herein. 


ANALYSIS 

Equations are derived that form the basis for the construction 
of charts describing the compressible air-flow process and hence the 
variation in pressure occurring in constant-area passages wherein 
simultaneous friction and heat-transfer effects are present. The 
charts are applicable to arbitrary heat-input distribution along the 
passage length. 

When the heat -flux distribution along the flow passage is not 
explicitly given but instead the passage dimensions, wall temperature, 
and air-flow parameters are given, a simple method of determining 
the temperature rise of the air as a function of the distance along 
the passage is desirable. From this method, the heat-input distri- 
bution for application in the method of pressure -drop determination 
described in the preceding paragraph may be obtained. Equations are 
derived for the temperature of the air as a function of distance 
along the passage wall for two cases : 

(a) Constant heat input along the passage 

(b) Constant wall temperature 

These equations, which are based on the recent heat-transfer corre- 
lation of references 5 and 6, are plotted to give convenient charts 
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for determining these air-temperature distributions as a function 
of passage length. 


Air-Pressure Drop with Arbitrary Heat -Flux 
Distribution along Passage Length 

The differential momentum equation for the one -dimensional 
steady- state motion of a compressible fluid in a constant-area pas- 
sage under the simultaneous influence of friction and heat-addition 
is 


d(mV + pA) + dDp = 0 (l) 

(All symbols are defined in appendix A.) In reference 7 it is shown, 
from energy and continuity considerations and the perfect gas law, 

that the total-momentum parameter ~ ■ + (hereinafter denoted by 

mA/gRT 

the symbol Mp) is uniquely related to each of the flow parameters 
M, p/P, pA/mf/gRT, t/T, and V/a/ gRT for any value of ratio of 
specific heats 7 corresponding to the total temperature T of 
the fluid. Eence, when the total-momentum parameter Mp and the 
total temperature T of the fluid are known, all the fluid-flow 
conditions (P, p, t, and V) are uniquely specified for a given 
mass flow per unit area. For convenience, the relations among the 
flow parameters are given in appendix B for the case of 7 equal to 
a constant. 

The variation of the total-momentum parameter Mp during the 
flow process involving friction and heat addition is due to the var- 
iations of both the total momentum mV + pA and the total tempera- 
ture T. Differentiation of Mp with respect to these two variables 
for constant mass flow gives 

dfcSjj&l 1 „« (2) 

mViRT * 1 

The differential drag force due to friction is 

„ 1 Tr 2 A 4 dx 

<% = F 2 pv A "D T 
6 


(3) 
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or, more conveniently, 

dDp = 2FmV (4) 

Combination of equations (l), (2), and (4) gives 

2_5L_ 

^2 = Fdx 1 dT f5 % 

Mp Mp D e " 2 T K 1 


The parameter cp is here introduced and defined by the equation 



VgRT 


( 6 ) 


The integration of equation ( 6 ) to give cp as a unique function of Mp 
and 7 is made in appendix C (equation (C3)). A plot of cp against 
Mp for air (7 = 1.40) is presented and discussed later in the section 

PRESSURE-DROP CHARTS. 


In terms of cp, equation (5) may be written 

Fdx \ 

_£e 1 dT 

dCp/ Mp 2 T 



(?) 


Integration of equation (7) between any two stations in a flow 
passage gives 


vM P,2 


"p,! 


Fdx\ 


D c 


)«£ _ 1 
/ 2 


dcp / Mp 


logp 


Ti 


( 8 ) 


where subscript 1 denotes flow conditions at an upstream station in a 
flow passage and subscript 2 , at a downstream station. 


The simplification can be made that over short distances of tube 
length Ax, cp changes linearly with x, that is, dx/dcp is closely 

given by — = 
b J Acp qp ? - cp 1 


where cpg corresponds to Mp ,2 an< ^ '■Pi 


cor- 


responds to Mp^q. With this simplification, the integration of equa- 
tion ( 8 ) may be formally carried out to give 
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F Ax 


\1 - 


cpi -cp 2 y ioge 


. Hz 

l°S e 


( 9 ) 


Plots of Mp^2 from equation (9) have been prepared for air (7 = 1.40) 
for increment in friction-distance parameter F Ax/D e = 0.04 and 0.01 
over a wide range of values of Mp,l and T 2 /T 1 . When the heat input 
per pound of air AH/mg in any increment of passage length Ax is 
known, the temperature T 2 may be obtained from air-enthalpy tables 
(reference 8 ); hence the plots of equation (9) permit calculation of 
Mp,2 from the known values of Tp, and AH/mg. The use of 


these plots for calculating pressure drop in heat- exchanger passages 

cLx Ax 

and the validity of the assumption — = - 

J dcp cp _cp 


are discussed later in 


2 T 1 

the sections PRESSURE-DROP CHARTS and Accuracy of Pressure-Drop Charts. 


For the special case of adiabatic flow (dT = 0), equations (5) 
and (6) give dcp = -Fdx/D e . A brief analysis of the adiabatic process 
with friction is presented in appendix C. 


Air-Temperature Variation along Flow Passage 
for Two Special Modes of Heat Input 

The air-temperature variation along the flow passage is herein 
related to the pertinent flow and heat -transfer parameters for 
(a) constant heat flux along the passage and (b) constant passage- 
wall temperature. 

The differential relation equating the heat transferred from wall 
to fluid to the heat absorbed by the fluid is 

mgc p dT » hs(T w -T) dx = dH (10) 

The heat- transfer coefficient h in equation (10) is given by 


hDg 

^w 


= 0.023 



( 11 ) 


This heat- transfer relation was obtained in the experimental investiga- 
tions of references 5 and 6 , which covered a wide range of air-flow 
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and heating conditions, including Reynolds numbers from 10,000 to 
500,000 and average wall temperatures up to 3000° R. 

In the analysis, the variation of the specific heat of air at 
constant pressure with temperature must be known. The empirical 
relation used is 

c p = KT 2 / 15 (12) 

which was found to closely approximate the actual variation of Cp with 
T ovei the range of temperature from 700° to 3000° R. 

From equations (ll) and (12) and the continuity relation, the 
following relation, of use in the analysis, is obtained: 


Ah 

mgc T 


= 0.023 


Ap, 

\jngD. 


0.2 


, T ») 


2/3 


. 0.6 


'p,w^wy 


(13) 


Constant heat flux. - The case of constant heat flux is character- 


ized by 


dE = ^t 
dx L 


(14) 


From equations (10) and (14) and the definition of D e , 


H, 


— mgCpT 


= 4 


Ah 

mgc T 


% > 

s. T J 


(15) 


De 


The total amount of heat added to the fluid in the over-all passage 
length L can be given as 


Ht - m S c p,av^ex " ^en) 


(16) 


or 


H+ 


®® c p , ex^ex 


C P> av fi _ 
c p,ex \ ^ex 

c P? av 
c 

p,ex 


( 17 ) 


With the use of equation (12), the ratio 
follows : 


may be evaluated as 
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L ex 


c p dT 


'P,av 

: P,ex 


L en 


15 i 


~p,ex v A ex 


< T ra - T en> * 17 Vt^/S 



From equations (17) and (18), 


17/15 


ex 


- T 


17/15^ 


en 


T - T 
ex en 


(18) 


H+ 


15 


“^^^ex^ex ^ L vex/ 


\17/l5 
1 en \ 


(19) 


Equations (13) and (15) are applicable at any cross section of the flow 
passage. When equations (13) and (15) are applied to the flow-passage 
exit, Cp and T are replaced by c p,ex and T ex , respectively, and 

the subscript w is replaced by the subscript w,ex. Hence, from 
equations (13) and (15), 


H + 


5 


’ °' 092 


0.2 


. 0 . 6 , 


"w, ex 


L ex 


\ c p,w,ex^w,ex y 


w.ex 


■J 


2/3 


L w,ex 

T 

ex 


( 20 ) 


(A\x 


From equations (19) and (20), 

\ 0.2 


w^exi 


mgD 


L_ 

D 


15 

17 


17/15 


1 - 


■ L en 

V. T eX/ 


0.092 


h L, 


ex 


\^pw,ex^w, ex/ 


f- s A-f 


T / V T 
en/ K± ' 


en\ 
w,ex/ 


2/3 


X w,ex 
^en / 


T 


en 


ex 


- 1 


( 21 ) 


Equation (21) relates .the over-all passage-distance parameter 

/ a \0.2 

( ^Pv ex 1 l 

\~lngD~J D~ " to ratio exit - "to entrance-air temperature 

T ex /T e n and the ratio of passage-exit-wall temperature to entrance-air 
temperature T v ,ex/Ten* For heat addition at constant heat flux, the 
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passage-exit-wall temperature is the maximum wall temperature. A plot 
of equation (21) is presented later for a constant value of 
Cp)i/k = 0.66. The total air- temperature rise in the passage AT 

corresponding to any given value of T e n and T w ,ex (maximum Tv 
for this mode of heat input) and the over-all passage-distance param- 
eter can be obtained from this plot. The total heat input cor- 

responding to Ten and AT can be obtained from air- enthalpy tables 
(reference 8). Inasmuch as for this case the heat input is uniform 
along the passage, the heat added in any distance along the passage 
is (x/L)H t ; from this value the air temperature at any cross section 
of the flow passage (corresponding to distance x) can be obtained from 
the air- enthalpy tables. 

Constant passage-wall temperature. - From equations (10) and (13) 
and the definition of D e for T w = K 



Integration of equation (22) gives 



(23) 


where, in the substitution of the integration limits, x is arbitrarily 
taken equal to zero at the passage cross section where T/T w equals 
0.20. This cross section may lie on a hypothetical forward extension 
of the given flow passage. The symbol x Q is the distance from this 
arbitrary zero point (where T/T w = 0.20) to the passage entrance, r 
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is the fraction of the actual passage length, and L is the actual pas- 
sage length. When r = 1, equation (23) gives the value of T/T w at 
the passage exit. A plot of equation (23) is presented later for 
c-pia/k = 0.66. From this plot of equation (23), the air- temperature 
variation, and hence the heat- input variation, along the passage length 
can be determined. 


PKESSUKE-DROP CHARTS 

Plots of equation (9) that apply to flow-passage segments for 
which the friction-distance parameter F Ax/D g = 0.04 and 0.01 are 
shown in figure 1. The total -momentum parameter at the exit of the 
flow-passage segment is plotted against the total-momentum 

parameter at the entrance to the flow-passage segment Mp,l for a 
range of values of total- temperature ratio Tg/'l across the seg- 

ment. In order to determine the variation in total-momentum param- 
eter (and hence all other flow parameters) with distance along a 
flow passage, the flow passage is divided suitably into segments of 
such length that the friction-distance parameter F Ax/D e is 0.01 
or 0.04 (or combinations thereof), and the total-momentum-parameter 
variation across each segment is obtained from figure 1. The dashed 
lines in figure 1 represent the condition where Mp^p = aud 

are drawn to facilitate use of the figure. The details of using 
figure 1 are illustrated later by means of examples I and III. 

When the FL/D e of the entire flow passage is not a multiple 
of 0.01 (and hence not a multiple of 0.04), at least one segment of 
the flow passage cannot be characterized by either F Ax/D e = 0.04 
or 0.01; hence for this one segment figure 1 is not applicable. In 
this case, the division of the entire flow passage into segments is 
so made that all of the segments except the first upstream one (the 
entrance station of which corresponds to the flow-passage entrance) 
are characterized by F Ax/D e = 0.01 or 0.04. The F Ax/D g of the 
first upstream segment will iieu.ce be less than 0.01. In this case, 
the approximation is made that the heat input in the first two 
upstream segments is added only in the second upstream segment and 
therefore the flow in the first upstream segment (of F Ax/D e less 
than 0.01) is adiabatic. The second upstream segment and all remain- 
ing segments are characterized by F Ax/D e = 0.04 or 0.01 and hence 
figure 1 is applicable for these segments. The foregoing approxi- 
mation as to location of heat addition near the entrance of a flow 
passage is justified cn the basis that: (l) At the entrance of heat- 

exchanger passages the comnressibility effects are relatively 
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small, in which case the exact location along the passage of heat- 
addition is not critical; (2) the tube length corresponding to 
F Ax/D e less than 0.01 will be a small percentage of the over-all 
tube length in practical heat exchangers. The procedure just 
described is illustrated in detail by means of Example III . 

In figure 2, the parameter cp is plotted against the total- 
momentum parameter Mp for air. It is shown in appendix C (equa- 
tion (C2)) that for the case of adiabatic flow with friction 

cpi~ cpg = £ - x . The use of figure 2 in conjunction with equation (C2) 
®e 

for determining the variation in total-momentum parameter along a flow 
passage for adiabatic flow with friction (no heat transfer) is illus- 
trated by means of Example II. 

The relations (for air) of the total-pressure parameter, the 
static-pressure parameter, and the Mach number to the total -momentum 
parameter are shown in figures 3, 4, and 5, respectively. The vari- 
ation of Mp and T in a flow passage can therefore be translated 

into the variation of total and static pressure and Mach number in 
the flow passage. 


AIR-TEMPERATURE-RISE CHARTS 


The relation of the air-temperature rise to the maximum passage - 
wall temperature (T v>ex in this case) and the over-all passage- 
distance parameter for constant heat flux along the passage length 
is given in figure 6. The ratio of the fluid-exit to fluid- entrance 
temperature T ex /T en is the ordinate, the over-all pas sage- distance 


parameter 


('A \ly eX \ 0 * 2 Xi 

, 1 — — is the abscissa, and the ratio of the maxi- 

ymgDg J D e 

mum passage-wall temperature to the entrance-air temperature 
T w ex/ T en is t ^ Le curve parameter. 


The ratio of the fluid to wall temperature is related to the 
distance parameter in figure 7 for the case of constant passage-wall 
temperature. In the construction of figure 7, the distance x was 
arbitrarily taken equal to zero when the ratio of air to wall tempera- 
ture equals 0.20; hence, the distance parameter in figure 7 must be 


^ w_\ 0 ’ 2 

\mgD e/ ' 


Xq +rl/ 
\ Dp / 


taken as 


where Xq is the distance of the passage 
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entrance from a hypothetical station in the forward extension of the 
passage characterized by T/T w = 0.20, L is the over-all passage 
length, and r is the fraction of the passage length at which the 
value of t/t w is to be determined. An illustration of the use of 
figure 7 in the solution of a pressure-drop problem is given in 
Example IV. 


EXAMPLES ILLUSTRATING USE OF PRESSURE-DROP CHARTS 

Example I 

The use of figure 1 is illustrated for flow with friction and 
heat addition when FL/D e is an exact multiple of 0.01. 

Given: 

(1) FL/D e = 0.270 

(2) mg/A =32.2 (lb/(sec)(sq ft)) 

(3) T en = 600 (°R) 

* 

(4) P en = 3000 (lb/ sq ft absolute) 

(5) The heat-input AE/a variation along the flow passage in 
terms of the friction- distance parameter Fx/D e is 


Fx/D e 

ah/a 

''"AH/mg 

2 T/T en 

0 

— 


1.000 

.04 

1018 

31.62 

1.218 

.08 

925 

28.72 

1.415 

.12 

848 

26.34 

1.594 

.16 

786 

24.40 

1.757 

.20 

713 

22.15 

1.904 

.24 

655 

20.35 

2.037 

.25 

157 

4.87 

2.069 

.26 

165 

5.13 

2.101 

.27 

157 

4.89 

2.133 


•*-The value of AH/mg is obtained 
by dividing AH/A by mg/A 
(item, 2) . 

^The value of T/T en is obtained 
from the air-enthalpy tables (ref- 
erence 8) for the given T e n 
(item 3) and enthalpy rise AH/mg. 
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Determine: The total and static pressures at the passage exit. 

The procedure is as follows: 

(6) From item (5) the values of F Ax/D e and T 2 /T 1 for the 
segments indicated in figure 8(a) are 


Segment 

F Ax/D e 

V T i 

a 

0.04 

1.2180 

b 

.04 

1.1617 

c 

.04 

1.1265 

d 

.04 

1.1023 

e 

.04 

1.0837 

f 

.04 

1.0699 

g 

.01 

1.0157 

h 

.01 

1.0155 

i 

.01 

1.0152 


(7) From items (2), (3), and (4), 

(■ — 5^— \ = 2.955 

•mr./gRT J en 

(8) From item (7) and figure 3, 

( ^ _ 3^126 = M,, -i (for segment a) 

Vm,/^T/ en ^ 

(9) From items (6) and (8) and figure 1(b), 

Mp,2 (for segment a) = 2.804 

This step is traced out in figure 9 by the line 1-a. 

Figure 9 is a reproduction of figure 1 and is inserted merely 
to aid in the explanation of the method of solution. 

(10) Mp^i (for segment b) = Mp ^ 2 (for segment a). 

Again from figure 1(a) and items (6) and (9) 

Mp,2 (for segment b) = 2.568 
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The operation of obtaining item (10) can be performed as illus- 
trated in figure 9. Proceed horizontally from point a to the 45° 
dashed line (point 2). The abscissa at point 2 is then Mp^i for 

segment b. Proceed vertically downward to the appropriate 

value for segment b (point b). The ordinate for point b is then 
Mp ^2 f° r segment b and also Mp ? i f° r segment c. 

(11) The operations indicated are repeated in step-by-step fashion 
for each of the segments a through i. For segments g, h, and i, 
use is made of figure 1(d) applicable for F Ax/D e = 0.01. These 

operations are traced in figure 9 where for each segment the point with 
ordinate and abscissa values that give Mp^g and Mp,i f° r se S~ 
ment is marked by the letter designation for the segment. (For 
example, point d in fig. 9(a) gives Mp,2 and M p,l for saS 111611 ^ &•) 

Although only the value for segment i is needed in this problem, a 
tabulation of the values of Mp^g is presented: 


Segment 

Mp,2 

a 

2.804 

b 

2.568 

c 

2.381 

d 

2.224 

e 

2.087 

f 

1.962 

g 

1.9313 

h 

1.8992 

i 

1.8670 


(12) For segment i, Mp^ 2 is Mp f or passage exit, as 

illustrated in figure 8(a). From figures (3) and (4) for Mp = 1.8670 
given in item (11), 



1.4865 


PA 


' m, s /gRT 


= 0.917 


ex 


(13) From item (5), T ex /T en = 2.133. From items (4), (7), and 
(12) therefore, 
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so that 



1.4865 

2.955 


A/2.133 


0.7349 


P ex = (0.7349) (3000) = 2205 (lb/sq ft absolute) 


(14) From item (12), 



0.917 

1.4865 


0.6167 


so that 


P ex = (0.6167) (2205) = 1360 (lb/sq ft absolute) 

(15) From item (12) and figure 5, the exit Mach number is 
M = 0.86. 


Example II 

The use of figure 2 is illustrated for adiabatic flow with 
friction. 

Given: 

(1) Mp = 2.500 at passage entrance (Mp, e n) 

(2) FL/D e = 0.300 

Determine: Mp at passage exit Mp^ ex 

(3) From item (l) and figure (2), cp at the passage entrance is 

cp en = 0.528 

(4) Equation (C2) of appendix C relates the change in cp between 
two stations in a constant-area flow passage to the friction-distance 
parameter between the two stations. Equation (C2) applied to the 
given problem can be written 

cp _ cp = HI: 

en ex p e 
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(5) From items (2) and (3), 

cp ex = 0.528 - 0.300 = 0.228 

(6) From item (5) and figure (2), 

Mp,ex = 2.182 


Example III 

The use of figures 1 and 2 is illustrated for flow with friction 
and heat addition where FL/D e is not an exact multiple of 0.01. 

Given: 

(1) FL/D e = 0.137 

(2) Wp,e n = 5 - 25 

(3) T en =. 800 (°R) 

The heat-input variation along the flow passage is specified and will 
he presented later in this example. 

Determine: Mp >ex 

The quantity FL/D e is not an integral multiple of 0.01. Hence, 
as illustrated in the sketch of figure 8(h), the division of the pas- 
sage into segments is so made that the one segment having a value of 
F Ax/D e other than 0.01 or 0.04 is the first upstream segment (seg- 
ment a in fig. 8(b)). Segment a has F Ax/D e = 0.007; segments h, 
c, and d have F Ax/D e = 0.04; and segment e has F Ax/D e = 0.01. 

The assumption is then made that the heat added in segments a 
and h is added only in segment h; hence the flow in segment a 
is assumed to be adiabatic, in which case figure 2 and equation (C2) 
are used to obtain Mp 2 for segment a, as illustrated in Exam- 
ple II. For the remaining segments, figure 1 is used, as illustrated 
in Example I . 
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(4) The actual heat- input variation along the passage and the 
heat input modified as described in the preceding paragraph are given 
in the following table: 


Segment 

F Ax/D e 

Actual 

AH/mg 

Modified 

AH/mg 

^Modified 

T 2/ T l 

a 

0.007 

9.76 

0 

1.000 

b 

.04 

51.93 

61.69 

1.313 

c 

.04 

80.41 

80.41 

1.300 

d 

.04 

101.60 

101.60 

1.280 

e 

.01 

33.30 

33.30 

1.070 


^The value of Tg/T^ across each section correspond- 
ing to the modified heat input is obtained from the 
value of T en (item (3)) and the modified enthalpy 

rise per pound of air AH/mg by means of the air- 
enthalpy tables of reference 8. 

(5) From item (2) (Mp >e n = Mp,l for segment a), equation (C2), 
and figure 2 for F Ax/D e = 0.007 (this step was illustrated in 
Example II ), 


Mp^g (for segment a) = 3.245 

(6) From item (5) \^,2 (f° r segment a) = Mp i (for seg- 

ment b)j and figure 1(b) for Tg/T-^ = 1.313 (item (4)), 

Mp^g (for segment b) = 2.803 

(7) For the remaining segments, the procedure is the same as 
given in item (6) for segment b (and as illustrated in Example I); 
hence only a tabulation of the results is presented as follows: 


Segment 

«P,2 

a 

3.245 

b 

2.803 

c 

2.425 

d 

2.100 

e 

2.017 
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(8) The value of Mp,ex for the flow passage is then given by 
Mp ^2 for segment e: 


The procedure for obtaining the remaining flow parameters at the 
flow-passage exit from the calculated value of Mp ; ex was illustrated 
in Example I. 


Determine: T/T en as a function of Fx/D e 

Pressure-drop measurements were taken in the experimental heat- 
transfer investigations of reference 6 and the data reduced to deter- 
mine the flow-friction factor F. The conventional correlation of F 
against bulk "Reynolds number was unsatisfactory. Improved correla- 


The best correlation line representing the data plotted in this 
manner is given by the equation 


This friction-factor equation is used in working out this example. 

P T 

(3) By taking — - = then from the friction- factor equation. 


Mp,ex = 2.017 


Example IV 

The use of recent friction-data correlation to relate Fx/D e 
and T/T en for constant passage-wall temperature is illustrated. 


Given: 


(1) T en /T w =0.320 


( 2 ) 



tion was obtained, however, from a plot of — against 

p w 




p av ^w 
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(4) With reference to figure 7, at the flow-passage entrance 
r = 0 and T/T w = 0.32 (item (l)). From figure 7, 



x 0 

% - 3 ' 67 


hence 3.67 represents the value of the distance parameter from a hypo- 
thetical station (in the forward extension of the flow passage) at which 
T/T w = the reference value of 0.20 to the entrance of the flow passage. 

(5) For the flow- passage exit, r = 1 so that from items (2) 
and (4) 



3.67 + 15.0 


18.67 


(6) From figure 7 and item (5), 


Toy 

= 0.777 
w 

(7) From items (1) and (6) with T qv 


1 + T 

en ^ X ex 

2 



0.549 


(8) From items (3) and (7), 

rx/B e ■ o.ow (^)°' 2 £ 

(9) The value of T/T w and Fx/D e at a cross section, say, 
nine-fifteenths of the distance from entrance to exit of the flow 
passage, is obtained as follows: 

r = 9/15 



rL 

D e 


= 9.0 
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( 

\mgDe, 

Hence, from figure 7, 


> 0.2 


Xq + rL 


D P 


= 3.67 + 9.0 = 12.67 


7T = 0-626 
1 w 


Inasmuch as x = rL, from item (8) the corresponding 
Fx/D e = 0.0285 X 9.0 = 0.257. 


(10) In a manner similar to that outlined in item (9), T/T v 
and Fx/D e can be determined at any passage cross section. The value 

t/t w t/t w 

of T/T en at the cross section is then given by 


A tabulation is made for various values of r 


'1A ' 

as follows: 


0.320* 


r 

/Apw \°‘ 2 x 
\ mgD e/ D e 

T 

T 

x v 

Fx 

D e 

T 

T 

^en 

0 

0 

0.320 

0 

1.000 

3/15 

3 

.426 

.086 

1.331 

6/15 

6 

.530 

.171 

1.656 

9/15 

9 

.626 

.257 

1.956 

12/15 

12 

.708 

.342 

2.212 

1 

15 

.111 

.427 

2.429 


A plot of Fx/D e against T/T en can be made to obtain T/T en for 
increments F Ax/D e of the passage equal to 0.04 or 0.01. The pres- 
sure drop in the passage is then obtained in the manner outlined in 
Example I or III. 


Accuracy of Pressure-Drop Charts 

In order to obtain an indication of the validity of the assump- 
tion used in obtaining equation (9) from equation (8) and hence of 
the accuracy of the pressure-drop charts, calculations were made for 
extreme conditions of heating and compressibility effects (l) by use 
of the charts and, (2) by numerical integration of the actual differ- 
ential equation (equation (8)). The conditions investigated and the 
results obtained are tabulated as follows: 
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M en 

^ex 

T ex/ T en 

Type of heat 
distribution 

Ap(exact) 
(lb/ sq ft) 

Ap ( chart ) 
(lb/ sq ft) 

0.30 

0.66 

2.04 

Constant wall 
temperature 

1111 

1107 

0.25 

0.46 

2.00 

Sine variation 
with distance 

241 

240 

0.30 

0.71 

2.00 

Sine variation 
with distance 

1273 

1274 

0.40 

0.85 

1.557 

Constant wall 
temperature 

1323 

1311 


For the first two cases in the preceding table, the charts for 
F Ax/D g = 0.04 were used to compute the chart pressure drop through- 
out the entire Mach number range from Me n to Me X . In the third 
and fourth cases, the chart for F Ax/D e = 0.04 was used up to Mach 
numbers of 0.67 and 0.65, respectively, and the chart for 
F Ax/D e = 0.01 was used for the remaining flow process. The agreement 
between the exact and chart results for any of these cases is within 
1 percent. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, March 30, 1950. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in the calculations and the 
figures: 


A 

C P 

D e 

% 

F 

S 

H 

s t 

h 

K 

k 

L 

M 

m 

P 

P 


cross-sectional area of flow passage, sq ft 
specific heat of fluid at constant pressure, Btu/(lb)(°E) 
equivalent diameter of flow passage, 4A/s, ft 
drag force due to friction, lb 


friction factor, 




ipV 

2 P D e 


mass conversion factor, 32.2 lb/slug 

heat added to fluid in distance x per unit time, 

Btu/sec 

heat added to fluid in total passage length L per unit 
time, Btu/sec 

heat- transfer coefficient between wall and fluid, Btu/(sec) 
(sq ft) (°B) 


constant 

thermal conductivity of fluid, Btu/(sec)(sq ft)(?R)(ft) 
total tube length, ft 
Mach number 

mass flow of fluid, slugs/sec 
total pressure, lb/sq ft absolute 
static pressure, lb/sq ft absolute 
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R gas constant for fluid (for air, 53.35), ft-lb/(lb) (°R) 

r fraction of length of given flov passage, 0 < r< 1 

s wetted perimeter, ft 

T total temperature of fluid, °R 

T w temperature of passage walls, °R 

t static temperature of fluid, °R 

V fluid velocity in flow passage, ft/sec 

x distance along flow passage as measured from passage entrance, 

ft 

Xq distance from hypothetical station characterized by 

T/T w = 0.20 to entrance of given flow passage, ft 

7 ratio of specific heats 

Ax length of flow-passage segment, ft 

p viscosity of fluid, lb/(ft)(sec) 

p mass density of fluid, slugs/cu ft 

Subscripts: 

av fluid conditions evaluated at average fluid temperature in 

flow passage 

en entrance of flow passage 

ex exit of flow passage 

w fluid conditions evaluated at average passage-wall tempera- 

ture T w 

1 entrance station of flow-passage segment 


2 


exit station of flow-passage segment 
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Parameters: 

% 

PA 

mfs/gRT 

PA 

mts/gRT 

V 

n/5rt 

CP 


total-momentum parameter 


static -pres sure parameter 


total -pressure parameter 


velocity parameter 


parameter, defined by dcp 


mV + pA 

m/s/gRT 


dMp 

V 


F Ax 


Fx 

D_ 


FL 

Do 



increment in friction-distance parameter 


friction- distance parameter 


over-all friction-distance parameter of given flow 
passage 


pas sage- distance parameter 


over-all passage-distance parameter of flow passage 
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APPENDIX B 


RELATIONS AMONG FLOW PARAMETERS 


For conven ienc e , the equations relating the parameters V//\/gRT, 
t/T, M, pA/mA/gRT^ and p/P to the total-momentum parameter Mp 

are herein derived. The assumption is made that 7 remains a 
constant. 


From the continuity equation and the perfect gas law, 

■ - oAV - H 

which can he written as 


PA \ 


V 


= Mp - 


V 


\ 


One form of the conservation of energy equation is 

2 


Y 


t = i . 2Ll_l ( , 

T ' \blm) 


From equations (B2) and (B3), 


(Bl) 


(B2) 


(B3) 


1 / V 'f y „ (_J_ _\ _2L 


N/gRT / 7 + 1 ' 7 + 1 


Mp 


= 0 


(B4) 


Thus 


V 


a/ gRT 


— is obtained as an explicit function of Mp 


V 


Vi RT 


- “p ± \!(yri “p 


2y 

7 + 1 


(B5) 
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From equations (B3), (B5), and the isentropic relation. 


7-1 



Then from equations (B5) to (B7), 



The parameter pA/m^/gRT may be obtained by direct substitution of 
equation (B5) into equation (B2). In the application of these equa- 
tions to a particular flow, the sign used in front of the radical should 
be plus if the flow is supersonic and minus if the flow is subsonic. 
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APPENDIX C 


METHOD OF INTEGRATION OF EQUATION (6) 

A brief analysis of the adiabatic compressible-flow process 
involving fric„tion only is herein presented together with the method 
of integration of equation (6). 


For the adiabatic flow of a compressible fluid with friction 
present, dT = 0 and hence from equations (5) and (6) 

• = -<§) 


(Cl) 


so that 


*1 





(C2) 


The variation of cp with Mp for either adiabatic flow or non- 
adiabatic flow may be obtained by substituting equation (B5), with the 
subsonic sign prefixing the radical, into equation (6) and integrating. 
The results following this procedure are 

r 


<P = i<Mp 


Mp 





where Mp and Mp^Q are the upper and lower limits of integration on 
Mpj <¥> and cpQ are the upper and lower limits of integration on cp, 
and has been taken equal to zero when Mp = Mp^ 0 . 


A large-scale plot of equation (C3) for y = 1.40 is presented in 
figure 2 for Mp^o = 1.8516 (the value of the momentum parameter at a 
Mach number of 1 with y = 1.40). 
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(a) 


M P,1 


Increment in friction-distance parameter (PAx/D e ), 0.04; range of total-momentum parameter at 

station 2, 1.92 to 2.52. 


Figure 1. - Charts describing compressible-flow process in a constant -area duct with simultaneous 
friction and arbitrary heat addition (y = 1«40)* (A 17- by 22-in. print of this fig. is attached.) 
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(b ) Increment in friction-distance parameter (FAx/De), 0.04; range of total-momentum parameter at 

station 2 f 2.52 to 3*01 • 


Figure 1. - Continued. Charts describing compressible-flow process in a constant-area duct with 
simultaneous friction and arbitrary heat addition (y = 1.40). (A 17- by 22-in. print of this fig. 
is attached.) 
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M P,1 

(c) Increment in friction-distance parameter (PAx/De), 0*04; range of total-momentum parameter at 

station 2, 3.01 to 3.48. 


Figure 1. - Continued. Charts describing compressible-flow process in a constant-area duct with 
simultaneous friction and arbitrary heat addition (y = 1.40). (a 17- by 22-in. print of this fig. 
is attached.) 
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(d) Increment in friction-distance parameter (FAx/De), 0.01; range of total-momentum parameter at station 2, 

1*852 to 2.109. 


Figure 1. - Continued. Charts describing compressible-flow process in a constant-area duct with simultaneous 
friction and arbitrary heat addition (y = 1.40). (a 17- by 22-in. print of this fig. is attached.) 
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( 0 ) Increment in friction-distance parameter (FAx/D e ), 0.01; range of total- 
momentum parameter at station 2, 2.105 to 2.226. 

Figure 1« - Concluded. Charts describing compressible-flow process in a constant- 
area duct with simultaneous friction and arbitrary heat addition (y = 1.40). 

(A 17- by 22-in. print of this fig. is attached.) 


1359 


NACA IN 2186 


33 



(a) Range of total-momentum parameter, 1.8516 
(choke) to 1.9488. 


Figure 2. - Relation between parameter cp and 
total-momentum parameter for air (y = 1.40). 
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Figure 2. - Continued. Relation between parameter cp 
and total-momentum parameter for air (y = 1.40). 
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Figure 2. - Continued. Relation between parameter qp 
and total-momentum parameter for air (y = 1.40) . 
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(d) Range of total-momentum parameter, 3*10 to 4.03. 


cp 


Figure 2. - Concluded. Relation between parameter cp 
and total-momentum parameter for air (y = 1.40). 
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(a) Range of total-manentian parameter, 1.852 (choke) to 2.418. 




Figure 3. - Relation between total-pressure parameter and total-momentum parameter 

for air (y * 1.40;. 
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(b) Range of total-momentum parameter, 2.41 to 3.83. 


Figure 3. - Continued. Relation between total-pressure parameter and total-momentum 

parameter for air (y ■ 1 . 40 ). 
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Figure 3* - Continued, delation between total -pressure parameter and total-momentum 

parameter for air (y ■ 1.40). 
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(d) Range of total-momentum parameter, 5*11 to 8.00. 

• - Concluded, Relation between total-pressure parameter and total -momentum 
parameter for air (y * 1.40). 


Figure 3 
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(a) Range of total -momentum parameter, 1.852 to 2.320. 


Figure 4 . - Relation between static-pressure parameter and total-momentum parameter for air 

(Y * 1.40). 
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(b ) Range of total-momentum parameter, 2.30 to 3.67. 

Figure 4 . - Continued. Relation between static-pressure parameter and total-momentum parameter 

for air (y = 1.40). 
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Figure 4. - Continued. Relation between static-pressure parameter and total- 
momentum parameter for air (y 3 1.40). 
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(d) Range of total -momentum parameter, 5*03 to 8.00. 


Figure 4 


Concluded. Relation between static-pressure parameter and total 
momentum parameter for air (y 3 1.40). 
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1.85 1.86 1.87 1.88 1.89 1.90 

M p 

(a) Range of total-momentum parameter, 1.8516 (choke) 

to 1.9500. 

Figure 5. - Relation between Mach number and total -momentum 
parameter for air (y = 1.40). 
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Figure 5. - Continued. Relation between Mach number 
and total-momentum parameter for air (y = 1.40). 
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(c ) Range of total-momentum parameter, 2.600 to 6.600. 


Figure b. - Concluded. Relation between Mach number 
and total-momentum parameter for air (y = 1.40). 
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Figure 6. - Variation of ratio of exit fluid temperature to entrance fluid temperature with over-all 
passage-distance parameter for various ratios of maximum (exit, in this case) wall temperature to 
entrance fluid temperature. Turbulent flow of air for constant rate of heat-input conditions. 

(A 17- by 22-in. print of this fig. is attached.) 
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Figure 7 . - Variation of ratio of fluid to wall temperature with passage- 
distance parameter for turbulent flow of air and constant passage-wall 
temperature conditions, 17. by 22-in. print of this fig. is attached.) 
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(a) Division of flow-passage length in Example I into segments a to i. 

Figure 8. - Illustration of method of dividing total flow-passage length into 
segments for use in pressure-drop charts. 
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(b) Division of flow-passage length in Example III into segments a to e. 

Figure 8. - Concluded. Illustration of method of dividing total flow-passage 
length into segments for use in pressu^ -drop charts. 
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(a) Increment in friction-distance parameter, 0.04. 

Figure 9. - Illustration of use of figure 1 describing compressible-flow process in constant- 
area duct with simultaneous friction and arbitrary heat addition (y = 1*40). 
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(b) Increment in friction-distance parameter, 0.01. 

Figure 9. - Concluded. Illustration of use of figure 1 describing compressible-flow process 

in constant-area duct with simultaneous friction and arbitrary heat addition (y = 1.40). g 
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(a) Increment in friction- distance parameter (FAx/De), 0.04; range of total- momentum parameter at station 2, 1.92 to 2.52. 
Figure I.- Charts describing compressible-flow process in a constant — area duct with simultaneous friction and arbitrary heat addition 
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Figure I.- Continued. Charts describing compressible-flow process in a constant - orea duct with simultaneous friction and arbitrory heat addition (ys|.4Q). 
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(c) Increment in friction- distance parameter (FAx/D e ), 0.04, range of total-momentum parameter at station 2, 3.01 to 3.48. 

Figure I.- Continued. Charts describing compressible - flow process in a constant - area duct with simultaneous friction and arbitrary heat addition (/ = l.40). 
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(d) Increment in friction-distance parameter (FAx/De)» O.OT, range of toto l-momentum parameter at station 2, 1.852 to 2.109. 

Figure I.- Continued. Charts describing compressible - flow process in a constant- area duct with simultaneous friction and arbitrary heat addition (y*l.40). 


NACA-Langley - 9-19-50 -1075 



1359 


230 


M 


P» 2 



(e) Increment In friction- distance parameter (FAx/D e ), 0.01; range of total- momentum parameter at station 2, 2.105 to 2.226. 

Figure I. - Concluded. Charts describing compressible - flow process in a constant- area duct with simultaneous friction and arbitrary heat addition (^=1.40). 
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Figure 6. -Variation of ratio of exit fluid temperature to entrance fluid temperature with over-all passage - distance parameter for various ratios of maximum (exit, in this case) wall temperature to entrance fluid 
NACA-Langicy - 9-19-50 - ms temperature. Turbulent flow of air for constant rate of heat-input conditions- 
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Figure 7.- Variation of ratio of fluid to wall temperature with passage-distance parameter for turbulent flow of air and constant passage-wall temperature conditions. 
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